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Abstract 

By implementing a genetic algorithm we search for stable vacua in Type IIB non- 
geometric flux compactiflcation on an isotropic torus with orientifold 3-planes. We 
find that the number of stable dS and AdS vacua are of the same order. Moreover we 
find that in all dS vacua the multi-field slow-roll inflationary conditions are fulfilled. 
Specifically we observe that inflation is driven by the axio-dilaton and the Kahler 
moduli. We also comment on the existence of one stable dS vacuum in the presence 
of exotic orientifolds. 
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1 Introduction 



Recently, there has been a huge interest in the search for classical de Sitter (dS) vacua within the context 
of superstring compactification. In the last few years, some constraints have been imposed by estimating 
the contributions to the effective scalar potential from each of the components that play a role in the 
compactification process. For instance, there is a no-go theorem, indicating the absence of dS vacua in 
compactifications threaded with standard NS-NS and R-R fluxes [1]. Further studies show, that by con- 
sidering a bigger set of allowed fluxes and more general structures for the internal geometry, it can be 
possible to find some stable dS vacua p]-ll5j. 

There are however some essential ingredients string compactifications should contain to enhance the 
chances of finding stable dS vacua: a negative curved internal manifold and a small number of moduli 
fields [Mil]. 

Here, we take the approach consisting in computing the scalar potential from a superpotential with all 
the above features. Specifically, a compactification on a negatively curved manifold with a superpotential 
W that depends at tree level on a small set of moduli is achieved by considering a type II string compact- 
ification on a six-dimensional isotropic torus in the presence of non-geometric fluxes |19H22j . 

These type of compactifications have become a fruitful region in the string landscape to look for stable 
dS vacua |23| I24j . For instance, by using an algebraic geometry approach |25fl30) . it has been shown that 
in the presence of non-geometric fluxes in Type IIB compactifications on /7j2 x Z2, there are stable dS 
vacua which are continuously connected to Minkowski non-supersymmetric vacua. Furthermore, recently 
the distribution of dS vacua and the presence of tachyon directions were studied in |3H |32] and a classifi- 
cation of stable vacua was shown in | 33j for a compactification on an anisotropic torus. 

As important as the search of dS vacua is, it is tempting to analyze these models in the context of 
inflation. Are there suitable conditions for slow-roll inflation in these models? Huge advances have been 
made in this field by considering general string compactifications |18| I33VI43) . One of the first models 
was studied in QJ, where it was proved that in the context of a standard vanilla compactification, the 
slow-roll cosmological parameters were much bigger than 1, excluding the possibility of inflation. Other 
models considered the presence of extra features, as more fluxes and extra stringy objects (orientifolds) 
yielding the construction of inflationary scenarios. However, ecently it was shown that single field inflation 
is forbidden [44] in the small-field regime. 

Explicitly computing the slow-roll parameters in specific models (as the non-geometric tori) requires 
numerical calculations. It is the purpose of this work to implement a genetic algorithm to search for stable 
dS vacuc[^and compute the corresponding slow-roll parameters, in order to establish if inflation is viable 
within these models. 

We find several stable vacua with positive and negative values at their minimum; interestingly we 
observe that the number of stable dS and Anti-de Sitter (AdS) vacua are of the same order. By com- 
puting the corresponding masses, we see that there is a hierarchy, pointing out the possibility to have 
multi-field inflation driven by a subset of the moduli. Specifically, we find that the inflaton is driven by 
the axio-dilaton and the Kahler parameter, with a direction approximately orthogonal to the s-goldstino. 
We observe that all stable dS vacua contain conditions for small multi-field driven slow-roll inflation. In 
particular we find that inflation is mainly driven by the internal volume and the real part of the dilaton field. 



^Previous advances towards this goal were shown in [45] . 
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We also find that supersymmetry is broken througli tlie complex structure at a trans-Planckian scale 
(with a high value of the gravitino mass for a realistic model) . Finally we show that a single field scenario 
for inflation is not allowed, in agreement with the results shown in |44| . We also comment on the existence 
of stable dS vacua in the presence of exotic orientifolds, an issue that has been so far avoided in literature. 
We believe that by considering exotic orientifolds it is possible to increase the moduli field space, where 
dS vacua can be found. 

It is important to remark that still it is not clear whether these models are constructed from a ten- 
dimensional scenario. The topic is still under study and great efforts have been made over the past few 
years |46| 147) . The outline of our paper is as follows: in Section 2 we fix our notation and study which 
moduli break SUSY while fulfilling the tadpole constraints and Bianchi identities. Section 3 is devoted to 
the implementation of the genetic algorithm to look for stable dS and AdS vacua. In section 4 we compute 
the slow-roll parameters and we show the possibility of inflationary trajectories for small fields. Also, we 
show that stable dS vacua are present in a background with exotic orientifolds. In section 5 we present 
our final comments and conclusions. 



2 Classical dS vacua 



We look for classical stable (A)dS vacua by considering a Type IIB compactification on a six-dimensional 
isotropic torus which is constructed by three identical copies of two-dimensional torus threaded 
with RR, NS-NS and non-geometric fluxes, under a Z2 symmetry given by the presence of 03-planes |21) . 
The effective = 1,D = 4 theory contains three complex moduli, namely the axio-dilaton S, a complex 
structure modulus r and a complexified Kahler modulus U. The contribution of the non-geometric fluxes 
to the superpotential is obtained by adding a term to the well known Gukov-Vafa-Witten |23j 



W{t,S,U)= I {F3-SH3-Q-U)An, 



(2.1) 



T6/Z2 

where Q ■ U = Q^^Umnrs]dx'^ A dx^ A dx^ is a 3-form and is the standard holomorphic 3-form. 

In such scenario the tree-level superpotential W depends on all moduli and after integration it is given 
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by 



(2.2) 



i=2 



where we have fixed the notation as follows: = (pi + iipi = {r, S,U} are the complex moduli fields and 
P($i) is a polynomial of cubic order on $1 with real coefficients Aij, 



(2.3) 



The real coefficients Aij in (2.3) are actually the integrated fluxes on the appropriate cycles of (see 
Appendix [B] for notation) given by 



(2.4) 
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The coefficients {ooi} correspond to the integrated RR-fluxes , {an} to NS-NS fluxes and {021} to 
non-geometric fluxes Q all over the 3-cycles on T^. Notice that for the isotropic case we are considering, 

3 

= 3a,i, for j = 0,1. (2.5) 

k=l 

We also consider the usual Kahler potential J(f of the form 

= -Y,Ailogi2^i), (2.6) 

i=l 

with positive Ai = {3, 1, 3}. Hence, if SUSY is broken by one or several moduli, the scalar potential reads 



(2.7) 



where the sum runs over the moduli that break SUSY and W = We 



ie 



We can see that it is possible to obtain positive values for 'f at its minimum. However, the existence 
of dS classical solutions depends on the presence of the non-geometric fluxes a2j- If they are turned off 
not only it is impossible to stabilize the Kahler modulus but also there is no possibility to find classical 
stable dS solutions since 

dip3 

which vanishes as ^/^s — )• 00 together with the whole potential, rendering the space-time to be flat. A 
similar situation holds in the absence of RR and NS-NS fluxes since in this case the scalar potential also 
contains a run-away direction, 

dy 1 

^ ~ 72' (2.9) 

that vanishes as '(/'2 — ^ 00, taking our model far away from the perturbative regime. Therefore, it is neces- 
sary to at least, turn on NS-NS, R-R and non-geometric fluxes, for the model to avoid the above run-away 
directions. 



The required conditions for dS or AdS vacua depends on which moduli breaks supersymmetry and how 
the vanishing of F-terms is compatible with Bianchi identities and tadpole constraints. 



2.1 Constraints 

As it is known, the presence of fluxes contributes with an amount of internal energy which must be 
cancelled by the addition of a negative tensioned object. These requirements lead to a set of constraints 
the fluxes must satisfy, as the R-R tadpole constraint ^ / A F3 + N = 16 and its T-dual counterpart 
F3 ■ Q = which after integration are given by 

3 

EE(-l)'«Oi«,V.) = 2(16-iV)5i„ fori = 1,2, (2.10) 

k i=0 

where the sum over k stands for the corresponding elements in the matrix Ojj, and is the number of 
exotic orientifold planes. An odd number of them are required at the fixed points for every odd flux 
supported in a cycle |18]; the extra constraints are encoded in the Bianchi identities. With the purpose 
to rewrite them, let us consider the vectors Aj^ := A^j^ given by 

^j,fc± - y^j{k)^"-j{k±l)-'^jik±l)^'^j{k±2))^ {'^■^^) 

(2.12) 
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where J = (j±, k) and M = (mib). The Bianchi identity Q ■ H-^ = decomposes as 

Af.Af: = {AfUAf),r^'' = ^, (2.13) 

for j 7^ y, diag{r]) = { — 1,1,1,1} and for the combinations M = (0+),J = {(2, 0+), (1, 3— )} and 
M = (3-), J = {(1, 2-), (2, 1+)}, while Q -Q = decomposes as, 

a ■ Cj = Ci^cjn,{-ir = 0, (2.14) 

for (i,j) = {(1,2), (1,3), (3,4), (2,4), with 

Cl = (a23, a22, 022) "2l)> 
C2 = («21) "^22; '^22) ^^23)) 

C3 = {a2o, all, all, 0^2), 

C4 = {al2,ali,ali,a2o)- (2.15) 
A formal derivation of the above constraints is found in 1211. 



2.2 SUSY breakdown 

Since we are not considering the presence of D-branes, there are no D-terms in the superpotential, and 
the breakdown of SUSY is achieved only through F-terms. Our model depends on three complex moduli 
$j, implying that there are 7 ways to break SUS^Vj^ It is therefore necessary to analyze all of them and 
check if they are compatible with the Tadpole and Bianchi identity constraints. 



A. Breaking SUSY through a single modulus. 

1. Consider that SUSY is broken only through the complex structure ^i, i.e., D^fiM^ = and DiW 7^ 0. 
The SUSY equations leads to the condition Im{P2Pz) = 0, which in principle can be compatible with 
the Tadpole and Bianchi identities, establishing a way to stabilize ^i. The solution of the SUSY 
equations of motion for <I>2,3 provides a stabilization for $2 but leaves unfixed the Kahler modulus 
$3. Its stabilization depends on the value of the scalar potential acquires at its minimum. We shall 
explore this possibility in more detail in the next section. 

2. Breaking SUSY only through the Kahler modulus <I>3 implies that the SUSY conditions -Di,2^ = 0, 
are fulfilled by any value of the moduli, even away from the minima. These two equations give the 
conditions 

ij^idiPi = 0, fori=l,3, (2.16) 

#l$25lP2 - 3'^2i^2 = 0. (2.17) 

Since we are assuming non- vanishing real moduli {(pi, ipi), the second condition implies P2 = 
(which in turn says that we have turned off all NS-NS fluxes. This gives a superpotential without a 
dependence on the axio-dilaton <I>2 and together with the first condition, it constraints the coefficients 
for -Pi, 2, 3 to be related by aoi = ^laij = \2a2j for each j and for any integers Ai,2- Then we obtain 
a violation of the Tadpole conditions. Since such cases lead to non-physical conditions, we exclude 
them from our analysis. 

3. Breaking SUSY through the dilaton modulus $2 leads to similar equations as the previous case and 
it will be also excluded from our analysis. 

®A complete solution of the supersymmetric cases is presented in |2U| with AdS vacua as expected. 
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Breaking SUSY through two moduli. 



1. Consider breaking SUSY through <I>2 and <I>3. Solutions to the equation DiW = imply that 
Pi($i) = -P2(*^*i) = -P3(*^*i)- The fluxes associated with similar power-terms in $i must have the 
same value. However this condition leads to aoj = aij = Xlfc '^2j J' ^^'^ violates the tadpole 
condition. 

2. Breaking SUSY through <I>i and <I>2 or through <I>i and $3 leads to Pj = for all i which also violates 
the tadpole condition. 

All these three cases are also excluded from our analysis. 

Breaking SUSY through all the moduli implies that all of them must be stabilized dynamically through 
reaching the minima of the scalar potential. The values of the corresponding fluxes must satisfy the 
Bianchi identities and the moduli must acquire positive values. Since there are no constraints to consider, 
this case must be taken into account in our search for stable dS and anti-De-Sitter vacua. 



In summary, we have found that there are only two ways to obtain stable vacua, compatible with 
Tadpole constraints, Bianchi identities and with the possibility to have a positive value at the minimum: 
breaking SUSY through the complex structure modulus <I>i or through all three complex moduli <I>j. It 
is important to remark that a generic flux configuration is not ad hoc compatible with these constraints. 
Therefore, we are interested in finding some specific flux configurations which fulfill such constraints. 



2.3 Inflation in non-geometric compactiflcations 

As we have seen, the existence of (stable) vacua with small positive or negative energy values depends on 
the way SUSY is broken. However, as shown by the expression for the scalar potential in (2.7), it also 
depends on the flux configuration we consider and on the vacuum expectation values for the moduli. 



Assuming the existence of such vacua (a task we face in the next section) it is of high interest to 
look for suitable conditions for the presence of inflation. A rigorous proof for that is quite complicated. 
However, we can infer some features about the existence of such conditions. 

Consider for instance a superpotential W = of the form 

3 

W = PY,Qi{^i). (2.18) 

for a constant complex P and for Qi being a polynomial of order n in for all i. The corresponding 
scalar potential can be written as 

e-^r = ^ {^\P?) + (2.19) 

i 

It is clear that the mass matrix is diagonal and the possibility to have all moduli masses at the same order 
is a viable option. Therefore, the possibility to identify a candidate for the inflaton within a subset of the 
moduli fields is excluded. Based on this argument, we think that a superpotential which does not look 
symmetric by interchanging the moduli fields $j (modulo the flux values aij) indicates the possibility of 
having inflation driven by a subset of the moduli . For instance, in KKLT-type models, it is customary 
assumed that dilaton and complex structure moduli fields are stabilized previous to the Kahler moduli, 
which acquire a vev by the presence of a non-perturbative term depending on them. In these models the in- 
flaton is identified with the Kahler and the axio-dilaton moduli or rather a linear combination of them |14) . 



6 



In our case, we have two (non-formally derived) reasons to believe that it is possible to identify a 
candidate for the inflaton with a single modulus complex field or with two complex fields . First, the 
superpotential W looks symmetric on $2 and <I>3 but not on $1, therefore we expect the inflaton to be 
identified with the complex structure modulus $1 or with a combination of $2 and $3. Second, we can 
show that in a generic scenario and with the expected conditions, the possibility that all masses are of the 
same order is very restricted, and in consequence, there is a big possibility that the model has a hierarchy 
on the masses. This points out the possibility to identify a candidate for the inflaton within a subset of 
the moduli fields. See Appendix [D] for details. 

In the next section we shall concentrate on two issues: 1) searching for stable classical dS and AdS 
vacua on which supersymmetry is broken and 2) studying the presence of suitable conditions for slow-roll 
inflation in the large and small field scenario. 



3 Searching for stable (Anti) De Sitter vacua: the genetic algorithm 

Let us start by looking for stable dS and Anti-de Sitter classical vacua. Notice that the existence of a 
minimum with all moduli stabilized depends on the values of the fluxes we are turning on. Looking for such 
vacua is a daunting task if one proceeds analytically. However it is possible to employ some computational 
algorithms which allows to find the existence of those vacua. 

We apply a genetic algorithrrj^ with a convergence criteria of 10^^*^. The numerical procedure is as 
follows. First a set of integer fluxes are randomly generated satisfying the Tadpole and Bianchi constraints. 
Then the genetic algorithm is employed in order to find the minima of the scalar potential. 

Searching for minima requires to fix a bound for all the fluxes [21] , which then acts as a natural 
cutoff, leaving the shape of the vacua distribution unaltered as the bound increases. However, a genetic 
algorithm does not allow to find metastable vacua since it searches for the lowest value of the potential 
in a neighborhood of a selected point. If another point is found in which the value of the scalar potential 
is lower than the previous one, the algorithm selects this point as the base to look for new minima. Once 
the searching fails to find a lesser value, the algorithm stops. Also, if a minima is found, the algorithm 
will not continue looking for other minima for the same potential. Therefore, we cannot be certain that a 
minima is metastable or global. 

In the next sections we shall employ the genetic algorithm to look for stable minima for the two 
cases derived in the previous section. For completeness, the algorithm code, with supersymmetry broken 
through all moduli, is presented in Appendix [C] 

3.1 Case 1: SUSY breakdown through the complex structure 

There are dS vacua solutions by breaking SUSY through the complex structure <I>i, as shown by the 
corresponding A-space distribution in Figure [T] It can be seen that most of the vacua are dS and only 
few of them are AdS. In both cases the energy vacuum values are much bigger than expected and most 

^ A genetic algorithm is a probabilistic optimization method to search global minima that mimic the metaphor of natural 
biological evolution. The computational procedure operates on a set of initial values for the moduli (population) which gives 
a potential solution by applying the principle of survival of the fittest to produce successively better approximations to a 
minima of the scalar potential. At each generation of the numerical procedure, a new set of moduli is created by the process 
of selecting individuals according to their level of fitness in the moduli space and reproducing them using feasible adaptation. 
This process leads to the evolution of moduli that are better suited to the minimum value of the scalar potential than the 
set of moduli from which they were created, so the algorithm successively improves the solution. In Appendix C we show 
the code used in this calculation. 
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Figure 1: Landscape for non- SUSY vacua with SUSY broken through the complex structure, (a) AdS, 
(b) dS. 



important, we observe that for all AdS vacua and almost all dS , the string coupling constant is very 
large, rendering our solutions unviable since they are far away from the perturbative regime where our 
approach is valid. Actually, there are only 3 cases for which gg is lesser than 1. However, the vacua 
distribution in terms of the internal volume (see Figure |2]) , lies in a region of small volume (small Kahler 
modulus) suggesting that the supergravity approximation is not valid precisely for those cases with Qs < 1. 
Therefore, for all found vacua in which SUSY is broken by the complex structure, there are not physical 
consistent conditions. 



3.2 Case 2: SUSY breakdown through all moduli 

Let us focus on a particular subset of the moduli space by considering only the case in which supersym- 
metry is broken by all moduli in the absence of exotic orientifold planes, i.e., with = 0. This forces us 
to consider only even fluxes. We shall consider the presence of exotic orientifolds at the end of the present 
section. 

After applying the genetic algorithm, we find two thousand vacua out of which only 15 have positive 
vacuum expectation values (vev's) and the rest are related to AdS vacua. For negative energy values, the 
landscape is shown in Figure |3^. However we find that almost all physically attractive cases are unstable. 
This follows from the presence of tachyonic directions in the eigenvalues of the mass matrix |49| [50] 

Mfj = Didjr + ^VoJ^ij, (3.1) 

where D is a covariant-Kahler derivative and Vq is the value of the scalar potential 'f at the minimum. 
The number of stable AdS vacua is of order 10 with half of them presenting vev's of order 10^^. In table 
[T]we show 5 representative cases. 



Concerning dS vacua, some of them correspond to solutions with a large string coupling constant gs where 
our approximations are not valid. The rest of them are related to large and small vev's. Focusing on those 
with small volumes, we report 10 different flux configurations leading to positive values of the correspond- 
ing classical scalar potential at its minimum, with small values of gs and small volume^ (see Figure [sja). 

^Notice that this does not mean that LARGE volume scenarios in presence of non-geometric fluxes are discarded. 
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Figure 2: Coupling constant versus volume of the internal space for all vacua. De Sitter vacua are marked 
as (8) while AdS vacua are marked as x. 



Table 1: Flux configuration for stable AdS vacua. 



Fluxes 


1 


2 


3 


4 


5 


o-oo 


2 


2 


2 


2 


2 


aoi 


16 


16 


18 


16 


20 


ao2 


8 


8 


2 


14 


18 


ao3 


4 


34 


20 


24 


38 


a23 


16 


16 


16 


16 


16 




2 


4 


2 


2 


4 


023 


16 


32 


18 


16 


40 



Flux configurations for stable dS vacua with the above mentioned properties are listed in Table [2] 



From these calculations we observe that: 

1. At least two non-geometric fluxes are required for the existence of stable (positive) vacua. We have 
not found a flux configuration with a single non-geometric flux leading to a stable —positive or 
negative valued— vacua. 

2. The order of stable dS and AdS vacua is the same (0(1)). As far as we know, this is the first case 
reported in which AdS is not favored against dS. 



Henceforth, within these 10 cases, all moduli are stabilized and their vacuum expectation values are shown 
in Table [sj together with the energy vacuum value which is of order 10^'' — 10^® in mass Planck units. 
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Figure 3: Landscape for SUSY breaking through all the moduli: (a) AdS vacua, (b) dS vacua 
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Figure 4: Coupling constant versus volume of the internal space. 
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Table 2: 


Flux 


configuration for stable dS 


vacua. 




Fluxes 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


aoo 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


aoi 


20 


16 


18 


18 


18 


20 


20 


20 


18 


16 


ao2 


40 


18 


20 


8 


34 


6 


8 


10 


30 


28 


ao3 


38 


16 


36 


20 


8 


8 


2 


36 


40 


2 


ai3 


16 


16 


16 


16 


16 


16 


16 


16 


16 


16 


"22 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


a23 


20 


16 


18 


18 


18 


20 


20 


20 


18 


16 



Table 3: Expectation value of the moduli at the minimum. 



Solution 








(V'2) 


(03) 


(^3) 






1 


0.01097 


0.0622 


36.728 


687.701 


1.2790 


854.393 


1.7331x10" 


-8 


2 


0.0219 


0.0661 


50.4234 


677.226 


0.2078 


861.222 


1.7830x10" 


-8 


3 


0.0132 


0.0530 


64.6363 


654.306 


1.4210 


856.003 


1.2476x10" 


'8 


4 


0.0064 


0.0702 


23.0519 


314.462 


1.1425 


416.45 


1.6714x10" 


-7 


5 


0.0049 


0.0735 


30.5774 


410.23 


7.1069 


512.971 


7.0674x10" 


'8 


6 


0.0078 


0.0561 


89.8675 


455.271 


0.1171 


555.596 


7.1009x10" 


-8 


7 


0.0075 


0.0639 


16.9536 


433.314 


0.5245 


501.969 


9.8853x10" 


-8 


8 


0.0020 


0.0700 


15.2937 


336.459 


0.9326 


432.311 


1.7225x10" 


-7 


9 


0.0285 


0.0658 


65.7282 


543.83 


0.9524 


673.968 


8.0765x10" 


-8 


10 


0.0106 


0.0691 


11.5051 


373.934 


0.0187 


499.693 


6.4109x10" 


-8 



The moduli masses are obtained by diagonalizing the mass matrix Mfj from Eq. (3.1 ) and their eigenvalues 



Mfj are given presented in Table |4| where is the squared mass of the eigenvector X constructed as a 



linear combination of the moduli fields (pi and ipi with X 
matrix P which diagonalizes the mass matrix by : 



PX and = {(</.i,V^i)}f^i 
P^^APP, is given by 



The 



P 



-0.8854 
-0.4648 
0(10-6) 

0(10-5) 
0(10-5) 
0(10-5) 



0.4648 
0.8854 

0(10 

0(10 

0(10 

0(10 



0(10-7) 
0(10-6) 
-0.1439 
0.0094 
0.4191 
0.8964 



0(10-5) 
0(10-5) 
0.0673 
0.2361 
0.8813 
-0.4036 



0(10-6) 
0(10-5) 
0.9855 
-0.0715 
0.0137 
0.1526 



0(10-5)- 
0(10-6) 
0.0577 
0.9690 
-0.2177 
0.1009 



(3.2) 



with a real redefined moduli field Xi given by 

Xi = PijXi- (3.3) 

The corresponding fields are approximately separated into two sets (cpi, ipi) and {(pi, ^i) with i = 2, 3. For 
i = 1, 2 we observe that j = 1, 2, while for i = 3, . . . ,6, the index j run over {3, . . . 6}. If we consider the 
fact that the masses for the complex structure moduli are much larger then clearly we have a hierarchy 
on masses. It seems that once we find stable dS vacua, all of them present such hierarchy. We shall 
extensively use this feature in the next section to find suitable conditions for inflation. 
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Table 4: Mass Eigenvalues. 



Solution 






■01 




)i>2 

<P2 




Mi 




]i>2 

03 








1 


6.354x10" 


- 


1.251x10" 


-0 


1.153x10" 


Q 

-o 


1.165x10" 


Q 

-o 


1.153x10" 


Q 

- o 


1.154x10" 




2 


1.110x10" 


-4 


1.908x10" 


a 
u 


1.179x10" 


c 

O 


1.191x10" 


a 
o 


1.179x10" 
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3.3 S-goldstino direction 

We are now ready to compute which modulu^ actually breaks SUSY: consider that for a particular set of 
vev's for all the moduli, we have a SUSY configuration, implying that all Kahler derivatives vanish for all 
the moduli. By departing from such values, SUSY is broken, rendering some of the F-terms to be different 
from zero. Since the F-terms are given by 

< F' >=< e^V-WJ^'~' >, (3.4) 

with the Kahler potential being diagonal, this only happens if < > becomes non-zero while the 
other are kept null, otherwise, as we have seen. Tadpole constraints and Bianchi identities are not fulfilled. 
Therefore, we conclude that the complex structure breaks SUSY. This is confirmed by computing the 
F-terms which for the complex field $2,3 are of order 10~^ while F^ is of order 10^. 

In the same line of reasoning, the gravitino mass 

m3/2 = (3.5) 

can be evaluated at the minimum of the scalar potential. Table [5] shows the values for the gravitino mass 
and the F^ term vev. 



Table 5: Gravitino mass and SUSY breaking scale. 

i 2 3 4 5 6 7 8 9 10 
TO3/2(xlO"'^) 1.796 2.143 1.701 3.217 2.282 2.049 2.449 2.293 3.017 3.028 
Fi-term 577.29 678.42 490.70 308.75 376.77 333.16 349.25 321.70 664.50 352.33 



We see that: 

1. The s-goldstino direction given by the F-terms, is approximately given by F^ indicating that it points 
towards the Kahler derivative with respect to <I>i ( i.e., the s-goldstino direction is identified with 
the complex structure direction). 

2. The gravitino mass and the scale at which SUSY is broken present very high values, around 10^^ 
GeV and 10^^ GeV, respectively. They are still far from being realistic. We briefly comment on this 
in section 5. 

^From now on, we shall refer to the redefined moduH fields $ as $. 
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Our next goal is to study if there are conditions for slow-roll inflation for small and/or large fields and 
determine if there is a single-field or a multi-field scenario for inflation. We also want to see if there is a 
preferential direction for inflation with respect the s-goldstino. 



4 Non-geometric slow-roll Inflation 

The presence of a hierarchy on the masses suggests a natural candidate for the inflaton, either composed 
by the complex structure ^'i or by a combination of the dilaton ^2 and the Kahler moduli $3. In order 
to establish which case is supported by our calculations we must consider small or large field scenarios. 
Since the last yields to a proper construction of a stable inflationary trajectory, which is quite complicated 
to find, we shall focus our analysis to small fields and leave some qualitative comments on the large-field 
regime after numerically computing the slow-roll parameters. 



In the small-field scenario, the complex structure $1 cannot be identified with the inflaton since its 
mass at the minima is much larger than the Hubble constant Vq. This leads to a possible multi- field infla- 



tion composed by <I>2 and $3. In such case, the complex structure would deca} ^ before the rest of moduli 
by which the goldstino would play the role of a Polonyi fielcj^ Notice as well that at this point there is 
still the possibility that inflation is driven by a single real field. We shall shortly see that this is not the case. 



4.1 Slow-roll inflation 

Despite the generalization to multi-fields it is still possible to have slow-roll inflation by defining an 
appropriate parameter that must be smaller than one. These are given by: 

= E^. = E(i^)"«i. (") 

More over, there must be conditions on which inflation ends, reflected as conditions upon which the slow- 
roll conditions are violated. 



In this case, the slow-roll parameters depend on 4 real fields and finding explicit and stable inflationary 
directions seems to be a difficult task. Instead we look for simple conditions for inflation to be present 
by projecting the slow-roll parameters in 2-dimensional planes where 2 real moduli are kept constant. If 
inflationary directions do not appear in such projections we could not conclude that multi-field inflation 
is absent. But if we are able to find some inflationary directions in these two-dimensional projections we 
can assure that at least we have a subset of all possible inflationary trajectories. Projections for the e 
parameter are shown in Figure [5] where the two of the moduli are taken constant at their stabilized values. 
Similar plots are obtained by keeping two of the moduli at constant values away from the minimum of the 
potential. 

From the above figures we conclude the following: 

1. The scalar potential is quite flat with respect to 03 and ip2 while the others are kept constant. 

^°Since we are not considering the presence of D-branes, the decay of the moduli fields via the interaction with other forces 
is absent in our model. The reader must understand this point as the possibility to decay if coupling terms in an effective 
action would be present. 

^^We do not have so far a solution for the cosmological moduli problem [511 152| . 
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e(>p2.(l'3) e(<!f.<l>2) 




200 400 600 SOO 50 100 150 200 



e) * f) 

Figure 5: Projections of parameter e in the plane a) V2-V'3 t)) 4^2-4'3 c) (p2-i^2 d) 02-'03 ^) 'p2-<Pz f) ^s-V's- 
The rest two real moduli are kept constant at their vev. The red circle denotes the position of the vev's 
for the two real fields. 



14 




2. There are important variations on the r] parameter as 02 and ip^ moves from the minima. 

3. It is possible to quahtatively notice the existence of inflationary trajectories on which inflation ends 
by facing zones where e reaches values of order 1. 

4. It seems that a combinations of (j)2 and ip^ would be a representative e-inflationary trajectory. 

5. Since V'a parametrizes the internal volume, we see that for an inflationary trajectory driven by 
variations of 'i/'s, far from the minima (still in the small- field scenario), the internal volume should 
be very small and it increases its size as it runs towards the minima of the potential. Therefore, 
inflation would be driven as the internal space increases its size. 

Let us now project the matrix terms rjij in two-dimensional planes and see if there are still inflationary 
directions in the field space. We select the matrix elements concerning only the fields (j)2 and ip^ since we 
have already see that e-inflationary trajectories exist. The projections are shown in Figure [6j from which 
we observe that there are conditions for inflation to exist and qualitatively trajectories for inflation to end. 



In summary, we have found that there are suitable stable inflationary directions. Particularly we have 
found that in the </>2-V'3 plane such directions exist. This leads us to the conclusion that at least for this 
case, the inflaton direction is (approximately) orthogonal to the s-goldstino. However, we should still 
show that inflation is sufficiently long to produce 60 e-folds. At this point this can be achieved either by a 
two-field inflation or by a single-field inflation, since we observe from Figure |5] that keeping constant one 
of the two real moduli (j)2 and -03 we still have inflationary trajectories. The number of e-folds required 
from phenomenological considerations to end inflation is given by: 



iV(0)= / ^^ ^60, (4.3) 
^J2ev{4>) 



4'end 
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a) 




Figure 7: a) A-i/^a vs {(j)2) b) A'i/'2 vs (tps). 



where i;^ is a linear combination of the moduh fields 02 and -03. In Figure [T] we show the plot of the 
e-folding values by keeping constant one of the two inflationary directions. As it is shown, it is impossible 
to cover the 60 e-foldings by small perturbations of the field, either by keeping constant (p2 or 03 . This is 
reflected from the fact that single- field trajectories joining the minimum with the region in which ~ 1, 
are sufficiently short to produce the required number of e-foldings. In other words, we find that single field 
inflation in the small- field scenario is discarded. This result is in agreement with the bounds computed in 

m- 

On the other hand, 60 e-folds can be achieved by many possible trajectories in the two-field inflation 
scenario, as shown in Figure |5] where a trajectory joining the minimum of the potential to a region in 
which e ~ 1 can be constructed by varying the two fields. It is still necessary to have a more detailed 
description about the explicit construction of such trajectories as well as a dynamical analysis about their 
stability. 

4.2 Stable vacua with exotic orientifolds 

Here we comment on the possibility to increase the number of stable vacua by allowing the presence of 
odd integer fluxes. This yields to consider exotic orientifolds in our setup [l8j. The presence of odd 
integer fluxes implies a modification on the Tadpole constraints. For the case of a single odd flux, the 
Bianchi identity is given hy ^ f A H3, = 15. But since all remaining fluxes are even, the identity cannot 
be fulfilled. For two odd fluxes (two exotic orientifolds), the Bianchi identity allows several solutions 
which are explored by searching the minima of the scalar potential with all the moduli stabilized. The 
corresponding landscape is shown in Figure [8] for which a further analysis reveals that only one case of the 
full-found landscape is actually stable. The moduli expectation values are given in Table [6] 



Table 6: Moduli and energy vev's. 




(02> (02) (03) (V'S) 




0.005 0.0251 


-134.55 5461.63 -70.085 1388.41 


2.8588x10-1° 



The corresponding non-vanishing flux configuration is given by the values: aoo = 2, oqi = 18, ao2 = 15, 
ao3 = 16, ai3 = 15, 022 ~ ^ '^23 = 18. Notice that odd fluxes are supported in different cycles as 
required and that this solutions implies the presence 62 03", one 03~^ and one 03 , which all together 
contribute to a total D3-brane charge of 15. 
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a) 



10'^ 



Figure 8: Landscape for N = 2 for a) AdS and b) dS. 

The moduli masses in the stable dS vacua are given in Table [7j from which we observe that a hierarchy 
is also present. Notice that masses for the U and S moduli are quite similar. This is a consequence of 
the mass term ^V^Jtij coming from the curvature of the space-time which is much larger than the second 
derivative DidjV. The gravitinos mass and the SUSy breaking scale are also quite high for realistic scenar- 
ios (around 10^^ and 10^^ GeV respectively). A similar analysis revels that slow-roll inflation conditions 
are also present for small fields, with the inflaton composed by a combination of the complex fields ^2 and 
^3. Since the features of this solution are quite similar to those without exotic orientifolds, we omit the 
detailed analysis. However we observe that considering the presence of odd fluxes increases the size of the 
field space related to stable dS vacua where slow-roll inflation is present. 



Table 7: Mass Eigenvalues for stable dS vacua for N = 2. 



Solution 




Ml 


Ml 

02 


Ml 






1 


1.341x10"^ 


1.445x10"'' 


1.906x10"^^ 


2.026x10"^'^ 


1.906x10"^'^ 


1.906x10"^^^ 



5 Discussion and Final Comments 

In this paper we have searched for stable de Sitter and Anti de-Sitter vacua in the context of Type IIB 
string compactification on a 6-dimensional isotropic torus threaded with NS-NS, R-R and non-geometric 
fluxes. These kind of scenarios contain all the essential ingredients for the presence of positive valued 
vacua. We have performed our searching by computing the scalar potential from the corresponding su- 
perpotential which depends on all moduli at tree level. This allows us to look for classical stable vacua 
within the context of a supersymmetric language. 

We have restricted our search to non-SUSY vacua and we have found that the single case compatible 
with Tadpole cancelation and Bianchi identities corresponds to the case in which SUSY is broken through 
all moduli, although the complex structure presents the higher F-teim. 

By implementing a genetic algorithm we have searched for stable vacua and we have found that the 
number of stable dS and AdS vacua are of the same order. With respect to dS vacua, we report 10 (similar) 
flux configurations with all constraints fulfilled in the presence of standard orientifold three-planes (i.e., all 
fluxes are even). In all of them, we find that the masses of the complex structure moduli are larger than 
the masses of the dilaton and of the Kahler moduli. We have also found that supersymmetry is broken 
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through the complex structure at a very high scale (around 10^^ GeV) with a gravitino mass of the order 
of 10^^ GeV while the vev's energy is of order 10^*^ GeV. We beheve that this issue could be improved once 
we incorporate interactions with other fields via the presence of D-branes. 

Concerning the existence of suitable conditions for inflation, we have found that in all the stable dS 
vacua we report, multi- field inflation driven by the axio-dilaton and the Kahler moduli is a viable scenario. 
Concretely, we have found that the inflaton is almost orthogonal to the s-goldstino direction and that the 
imaginary part of the dilaton and the internal volume are the real fields which contribute more significantly 
to the slow-roll parameter values. In particular we show that single field inflation is discarded in agreement 
with the bounds computed in fl^ . An inflationary direction driven by the above two real fields is a viable 
direction according to our calculations. However, specific directions requires more detailed considerations 
which are beyond the scope of this work. Finally we increase our field space by turning on odd fluxes. As 
it is well known this implies the presence of exotic orientifolds. In this case, we report one single case for 
a stable dS vacuum with the same characteristics as the vacua found in the absence of exotic orientifolds: 
a high scale for SUSY breaking, high gravitino mass, a hierarchy on the moduli masses which separates 
the complex structure from the axio-dilaton and Kahler moduli. Single field inflation in the small regime 
is also discarded. 

It is worth mentioning some issues which limited our study. First of all, the use of genetic algorithms 
makes impossible the search of Minkowski vacua since it works on a numerical approximation. Second, 
this method also excludes the possibility to conclude either a minimum is meta-stable or global. However 
we can assure that in case of having a lower minimum, it must be more than a 100 mass planck units 
away from the reported one. Metastable vacua close to the studied minimum are not considered by the 
algorithm. Finally, since the complex structure moduli decay previous to the inflaton field, it plays the 
role of a Polonyi field. So far we have not a solution for the cosmological moduli problem derived from 
this scenario. We leave this feature for further work. Finally, we should say that inflation driven by all 
moduli or by a single one is not discarded in a large field scenario. 
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A Scalar potential 



In Figure (|9| we show the plots concerning the minima of the potential for Case 1 with a positive valued 
minima in the absence of exotic orientifold planes (A^ = ) and for the single case reported in the presence 
of two exotic orientifold planes {N = 2). Observe that in both cases, the complex structure directions of 
the potential are not as flat as along other directions. 




Figure 9: a) Complex Structure, c) Dilaton, e) Kahler modulus for A'^ = 0, b) Complex Structure d) 
Dilaton f) Kahler modulus for N = 2. 

B Notation 

In this section we present for completeness the notation we follow in this work and its relation with the 
one commonly used in the literature. As usual the holomorphic 3-form can be written as; 

^ = (dx^ + Tdy^)A{dx^ + Tdy'^)A{dx^ + rdy^) (B.l) 
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2 

where J is the usual Kahler (l,l)-form. 



U = C4 + \e-'^ J = ^\JiZj\ 



(B 



where = x*^^. The orientifold 03 involution acts as (T(a;*) = — for i = 1, ...,6. The 3-forni with 
leg on each torus are; 

ao = dx^ Adx'^ Adx^ /3o = dy^ Ady'^ Ady^ 
ai = dx^Ady'^Ady^ /3i = dy^Adx'^Adx^ 
"2 = dy^Adx'^Ady^ = dx^Ady'^Adx^ 
03 = dy^Ady'^Adx^ P3 = dx^ Adx"^ Ady^ , 

with the normalization Jj,g aj A Pj = 5jj. So, the holomorphic 3-form is rewritten as 

3 

Q = ao + T^l3o + ^{Tl3i + T^ai), (B 

i=l 

and the closed 2-forms and their 4-forms duals are 

Ui = -dx^Ady^ Uk = dx^ Ady^ Adx^ Ady^ , (B 
where A u)j = 5ij. The NS-NS flux is written as; 

3 

-^3 = aio/?o + 01300 + ^ (aiiOj + ai2/3i), (B 

1=0 

while the RR 3-form flux reads 

3 

-^3 = aoo/3o + 003^0 + ^ (ooiai + ^02 A)- (B 

1=0 

Finally, the Kahler moduli can be written as 



fB 



The integrated non-geometric fluxes Q are given by the numbers a2j as shown in Table [8j 



Table 8: Fluxes in the duality invariant superpotential. 
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flux 


1 
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ooo 




$1 


F- ■ 


aoi 




$2 




ao2 




$3 


FafB'y 


ao3 
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Hijk 


aio 




$3 


Qf 


020 




$2^>i 




On 








O21JO21, 


oil 


$2^>f 




ai2 






Q-f iQ~p iQk 


022)022' 


"22 


^2^1 




Ol3 




$3^-? 




O23 
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C Algorithm Code 



clc 

clear all 

global aO al a2 a3 bO bl b2 b3 cO cl2 cl3 c22 c23 c3; 

max=1000; % Maximum number of vacua in the landscape 

Nu=10; % Cutoff for the fluxes 

nvar=6; % Number of variables in the scalar potential 

time_ga_parallel=0 ; 

sol=zeros (max , 7) ; 

f luxes=zeros (max, 14) ; 

fprintf('\n Minimization of moduli begins\n\n' ) ; 
for i=l:max; 
clc ; 

numero=randi (10) ; 

[aO al a2 a3 bO bl b2 b3 cO cl2 cl3 c22 c23 c3]=solucion(numero,Nu) ; 

lb=[0, 0,0, 0,0,0] ; 

Vl=a(x) \l. Function to minimize 

options = gaoptimset( 'Generations' ,4000, 'TolFun' ,le-20, 'PopulationSize' ,100, 'TolCon' ,le-4, . . . 

'MutationFcn' ,Smutationadaptf easible, 'StallGenLimit ' ,35, 'Display' , 'off , ' InitialPopulation' , initial, 'HybridFcn' ,@fmincon) ; 
startTime = tic; 

[x,fval,flag]=ga(Vl,nvar, [],□,□, □ ,1b, [] , [], options) ; 
time_ga_parallel = toe (startTime) ; 

sol(i,:) = [x, f val] ; 

f luxes(i, :) = [a0 al a2 a3 bO bl b2 bS cO cl2 cl3 c22 c23 c3] ; 

fprintf ('Iteration \y.g takes \y,g seconds. ( \y,g, \y.g, V/.g , \y.g , \y.g , \y.g , \7.g) \n' , i ,time_ga_parallel , sold , 1) ,sol(i,2) , . 
sol (i, 3) , sold, 4) ,sol(i,5) ,sol(i,6) ,sol(i,7)) ; 

end 

f igure(l) ; 

loglog(l./sol(:,4),-sol(:,7),' x '); 
xlabeK 'g_s ' ) ; 
ylabel('ln(\Lambda )'); 



D A window for inflationary conditions 

The existence of stable De Sitter vacua, as we have seen in the body of this manuscript, depends on the 
way SUSY is broken and the specific flux configuration we are considering. However, from the point of 
view of cosmology, it is interesting to look for extra conditions, namely: 1) Small values for the vacuum 
energy (positive or negative) and 2) conditions for the presence of inflation (slow-roll parameters for the 
effective scalar potential). 

Here we want to show that there are generic issues which allow us to expect that there are big chances 
to find specific models which allows the presence of inflation by assuming small values for vacuum energy 
(although this requirement is not necessary). 

Our strategy is to consider a subset of solutions on which the possibility to have an inflation driven 
by a subset of the moduli space is much bigger than its absence. This subset of solutions represents the 
"window" we are going to focus on. The idea is the following: we shall assume that it is possible to 
identify a subset of the moduh fields as the inflaton. This would require that the masses of such fields 
are larger than the rest. If we can show that this is the most probable scenario we will conclude that the 
identification of those fields with the inflaton is also a very likely situation. To do that, we will proceed the 
other way around, i.e., we will show that assuming all masses to be of the same order is a very restrictive 
option, while having the contrary represents the big subset of solutions. This argument does not imply 
the existence of inflation conditions, but open up the possibility to find them. 
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D.l Mass hierarchies 



We shall show that there exists conditions for inflation by probing that the absence of inflation is very 
restrictive. Our statement is as follows: If all masses computed from the scalar potential are of the same 
order, there would not be a hierarchy on the masses which otherwise would lead to the presence of a prefer 
inflationary direction driven by a linear combination of moduli fields. 

Therefore, if constraining the scalar potential terms such that all masses are of the same order in 
Planck units is not consistent, it would mean that all moduli fields would decay at the same rate and 
there would not be an inflationary direction on the potential, although a multi- field inflation driven by all 
moduli could be favored. 

Let us therefore assume that all masses computed from the scalar potential are of the same order. This 
statement restricts the terms of the polynomials Pj ($i). Since the square magnitude of these polynomials 
appear in the mass terms, their components must be of the same order since all of them depend only on $i. 

We shall concentrate on the diagonal terms in mass matrix (since this is enough for our purposes). Let 
us consider first the masses related to the moduli fields $2 and $3. It is easy to check that 



Mis,, = ^ ~ 2 M\ + 1 (D-1) 



r2 _5^r^2f^ -^---2 

for i = 2,3. If all terms in the polynomial contribute significantly to the mass, this means that 

0{< ajo >) ~ 0{< a^i^i >) ~ 0{< aj2^j >) ~ 0{< aj^^f >), (D.2) 

where < $1 > is the value of <I>i at the minimum of the potential (if exists). Consider a generic term 
djm^T- If such term is present it must be of the same order that o.j(m+i)^T'^^ ■ This implies that 

< ajm >~< «j(m+l)^l > • (D.3) 

Let us say that < $1 >~ lO^Mp, for a real T < 1 (since we are assuming small values on the scalar 
potential). Let us also assume that the real coefficients related to the integrated fluxes, are all of 

order 10''^ in mass Planck units. For the above two terms of the polynomial to be of the same order we 
find that 

0{< aj^ >) ~ 10^+^, (D.4) 

in Plank units. Therefore we would require a fractional amount of fluxes violating Dirac quantization. 
Hence, Pj($i) is restricted to be of the form 

Pii^i) = a^m^T. (D-5) 

for some m G {0, 1, 2, 3}. Notice that at least one of the polynomials must depend on $1, i.e., for some i, 
m must be different than 0. It follows then that both polynomials P2 and P3 must be of the same order 



at the minimum of the potential. From expression (2.3), we observe that this happens if 7712 = m^. 



Now, consider the mass term M^^^^, which assuming that 0(^2) ~ ^'(Ps) is given by 

^^1^1 ~ (2^ + l)l^l' + 2"^2 j^{i^l\P2f + i'l\P3\^ + i^2 Im {Wp2)+ti Im [WP^)). (D.6) 

This mass term would be of the same order as M^^^^^ and M^^^^^ if C'(V'2-P2) ~ ^(V's-Ps) ~ ^{Pi)- O^^^ 
possibility consists in taking a vanishing Pi, but this would imply the absence of NS-NS fluxes violating 
the tadpole condition. Therefore Pi must be different from zero, leading to the constraint that 

0{< $f ^ >) ~ 0{< V'a^r >)• (D-7) 
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Assuming that the order of < ■02 > is < $i with an arbitrary s, we are forced to say that mi = s + m2. 
However, since the numbers mi are the degrees on the polynomial coefficients, they must be integer num- 
bers. 

On the other hand, we also see from the tadpole condition that m\ — m2 = 1, 3. Therefore, a scenario 
where all masses are of the same order yields to a vacuum expectation value for the real part of the dilaton 
and the Kalher structure to be of the same order or 3 order of magnitude bigger than the vev of the 
complex structure. In summary, we find that in order to have all masses of the same order we must fulfill 
the following conditions: 

1. Small values for < $i >, which lead to small values for the scalar potential (since it is constructed 
by positive defined terms). 

2. All fiuxes of the same order. 

3. 0{< $2 >) ~ 0(< $3 >) ~ {0{< $1 >))^'^/^ 

This scenario in which inflation guided by a subset of all moduli is not present seems to be quite 
restrictive due to the condition (3). The possibility that inflation is indeed driven by a subset of the 
moduli is therefore the biggest one although we can relax the flrst two points. Relaxing point (1) would 
lead us to big values for the minima of the scalar potential. The second point can nevertheless be relaxed, 
but no so far from small values for the fluxes, otherwise we could depart from smooth geometries and 
backreactions must be considered. Notice that the mass terms above considered corresponds to elements 
of a non-diagonal mass matrix. Once it is diagonalized by computing the proper values, the order of the 
diagonal terms keep close to each other. 
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